Perturbative behaviour of a vortex in a trapped Bose-Einstein condensate 



Lyndon Koens 1 and Andrew M. Martin 1 

1 School of Physics, The University of Melbourne, Parkville, 3010, Australia 

(Dated: March 1, 2012) 

We derive a set of equations that describe the shape and behaviour of a single perturbed vortex 
line in a Bose-Einstein condensate. Through the use of a matched asymptotic expansion and a 
unique coordinate transform a relation for a vortex's velocity, anywhere along the line, is found in 
terms of the trapping, rotation, and distortion of the line at that location. This relation is then 
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I. INTRODUCTION 

Experimental realisations of vortices in Bose-Einstein 
condensates (BECs) [JHH] has spurred lots of research 
into the behaviour and effects vortices have on a BEC 
[5J [7] . Similar to the vortices Feynman predicted in su- 
perfiuid He II jS], vortices in BECs have quantised circu- 
lation, ensuring that the condensates phase is a single val- 
ued function (see Ref. [9] for details). However unlike in 
superfhud He II, the vortices in a BEC are sufficiently de- 
scribed in the mean field regime by the Gross-Pitaevskii 
equation [TUHH] with rotation [Eq. Q]. 
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where 'J is the condensate wavefunction, V tr (r) is the 
external trapping potential, l~i is the condensate rota- 
tion vector, g = Aixaf? /M is the interparticle interaction 
strength, /Lt(O) is the chemical potential, M is the particle 
mass, and a is the s-wave scattering length. 

In a stationary condensate it is energetically un- 
favourable to contain vortices, however when > f2 c = 
5Mn {Rj_/r c ) /2MR 2 ± (R x being the radius of the BEC 
perpendicular to the rotation vector and r c being the 
vortex core radius) a centred vortex becomes energeti- 
cally favourable [12]. Upon investigating the dissipation 
of such a vortex, the vortex's dissipation time was shown 
to depend on the condensate's temperature [T3]. 

In a trapped BEC a quantised vortex is rarely station- 
ary, usually moving or contorting. An off centred straight 
line vortex is known to precess at <j) — 3f2 c /5(l — pg/i?^) 
(po being the cylindrical radial coordinate of the vortex 
line) [12 114] . This straight structure is representative 
of a vortex in a pancake shape condensate, R±/R z > 1, 
where R z the condensate radius parallel to the rotation 
vector. For a cigar shaped trap, R±/R z < 1, however 
the vortex line bends. This result was shown numerically 
in Refs. [HI [16], seen experimentally by Rosenbusch et 
al. [T7], and reasoned physically by M. Modugno et al. 



|18) . A quantised vortex also supports wave perturba- 
tions, coupling with and modifying the normal modes of 
the condensate [HI [2U] or supporting helical wave struc- 
tures, like Kelvin waves [21], along its length p~0| [22] . 

With such a wide range of behaviour, a single set of 
equations that completely describes the vortex has been 
desired. Early attempts to derive this equation relied on 
the method of matched asymptotic expansion [231 124] . 
This procedure was extended by A. Svidzinsky and A. 
Fetter [25] to produce a set of equations that should gen- 
erally described all small perturbations from a straight 
line vortex. This derivation used a coordinate transform 
to account for vortex line bending, however in its exe- 
cution, these coordinates were not formally defined, and 
the outer and inner solutions matched through a pseudo 
'vector potential'. This gave rise to a modified set of 
equations to those that describe the vortex's motion and 
structure. 

This paper re-derives this procedure to determine the 
behaviour and structure of a slightly perturbed straight 
line vortex, using coordinate transformations and the 
method of match asymptotic expansion. Section [IT] justi- 
fies the relevance and procedure of such an expansion, 
and Section |III| defines the unique coordinate system 
used. An inner and outer solution are then found (Sec- 
IV and |v| respectively) and matched (Section VI) 



tions 



to give a relation between the perturbed vortex's motion 
and shape. This relation contains an unknown constant 
that is determined by comparing to known physical sce- 
narios (Section |VII[ ). This then allows simple results from 
the relation to be calculated and compared to numerical 



simulations (Section VIII). Finally in Section IX a gen- 
eral comparison between this work and Ref. [25] is per- 
formed. Illustrating the differences in methodology and 
general results. 



II. THE IDEA 

BECs with vortices within them have two natural 
length scales: the condensate length scale R± and the 
vortex core radius r c , with r c < R±. These two scales 
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suggest that a matched asymptotic expansion can 
be performed to determine the behaviour near the vor- 
tex core, length scale r c , far from the vortex core, length 
scale R±, and matched to give the systems full behaviour. 
Such an asymptotic expansion has been previously pre- 
formed for vortices in simplified scenarios in Refs. [23] 
and [21]. 

In the local coordinates of the vortex line the inner 
solution depends on properties of the vortex line, which 
when matched to the outer solution creates a relation 
between the shape and velocity of the line. This rela- 
tionship then gives a set of differential equations that 
describes the vortex line's structure and motion through 
space. 



III. THE LOCAL COORDINATES OF THE 
VORTEX LINE 



In a BEC the localized vorticity from a vortex, cre- 
ates a single vortex line. Each point on this line has a 
position and a velocity; therefore, the vortex line, at a 
given time, is described by two parametrized functions: 
a curve describing the lines shape and position, c(T), and 
a velocity vector, V(T). This parametrization allows for 
the transformation into the local coordinate system; the 
curve, c(T), determining the axes and V(T) representing 
the line's motion at a given point. 

For any curve c(T)], a tangent vector [Eq. j2j], a nor- 
mal vector [Eq. (3)], and a binormal vector [Eq. Q] are 
defined, see Fig. T These three orthogonal vectors form 
the basis to the coordinate system of the curve, known 
as the Frenet-Serret coordinates [2*7] . 



t = 



c'(T) 
|c'(T)|' 

c'{T) x (c"(T) x c^T)) 
|c'(T)| |c"(T) x c'(T)| ■ 

t x n. 
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At a given T, such a basis is equivalent to the basis of 
a Cartesian coordinate system. Allowing the respective 
t, n, and b vectors to be locally treated as z, x, and y. 

This local Cartesian behaviour extends to the deriva- 
tives in n and b directions. Translations in these direc- 
tions do not change the basis vectors, and therefore the 
local derivatives in n and b are identical to their Carte- 
sian counterparts (d n = d x , d% = Q\ , etc.). 

The curve progresses in t causing the local derivatives 
in t to not match its Cartesian equivalent. Using the 
chain rule, this deviation can be quantified, showing the 
first order derivatives in t to be 



d_ 
df 




Figure 1: Depiction of the local coordinate vectors t, n and b 
ofacurvec(T) = {T 2 /2, T 2 /2, T} (blue), with the coordinate 
vectors plotted for T = (red) and T = 1 (green). 



from which higher order derivatives can be constructed. 
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where ® represents the outer product and : represents a 
double contraction. 

In this case the structure of the vortex line is unknown. 
Therefore in order to evaluate these derivatives a pseudo- 
parametrization, that describes the local behaviour and 
can take any shape, needs to be used. 

Conveniently, any unique curve is defined through 
two terms: curvature [Eq. Q] and torsion [Eq. Q]. 
These terms describe how a curve bends and distorts 
through space, irrespective of its location and choice 
of parametrization. Curvature indicates how the curve 
bends in the t, n plane (with 1/k being the radius of 
a circle at a given point) and torsion indicates how the 
curve twists out of the t, n plane (how n and b rotates 
as T progresses): 



|c'(T)xc"(T)| 

|c'(T)| 3 
(c 7 (T) x c"(T)) ■ c"'(T) 

|c'(T)xc"(T)| 2 



(7) 
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Generally, a curve's torsion and curvature vary along 
its length, allowing for all possible 3D curve structures 
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to form. However, locally around any point the torsion 
and curvature is effectively constant, prompting the use 
of a modified helix for the pseudo-parametrization of the 
vortex line: 



c(T) = {a cos(T),a sin(T),6 T} . 



(9) 



Assuming the excitations on a vortex line are small 
compared with the vortex's overall structure the curves 
tangent vector will almost align with z. This condition 
can be enforced by adding appropriate linear terms to 
the parametrization, to make t|T=o ~ z - This condition 
also implies that the length of the curve is predominantly 
the length travelled in z (b is large); this simplifies the 
re-normalization of the parametrization, from T into the 
arc length of the curve s f=a z = bT. This is an easier 
parametrization to work with as curves parametrized by 
arc length have |c'(s)| = 1. 

Furthermore, helices have no variation in p. Hence 
for the complete curve to take any form, a constant 
pseudo-curvature k, between p and z, needs to be added 
to the pseudo -parametrization [k being defined as ps 
a/p"(s) 2 + z"(s) 2 ]. Keeping z(s)'s form [z(s) = s], p"{s) 
equals —As. Hence the pseudo-parametrization becomes 
[Fig- 1 



c(s) 



as, 



(10) 



where a and (3 are the small angles of deviation from z in 
the x and y directions respectively, i.e. t| = {a, (3, 1}. 

This parametrization makes the curvature, torsion, 
and local derivatives in s, to first order in a and (3, to be 
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This makes the gradient and Laplacian operators (V and 



'-) 



V = V 



d_ 

dx 



4 

dy 



z. 



(15) 



and 



V' 2 = V 2 - K 



dx 
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dx dz 
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where V has the standard Cartesian definition (V 
{d x ,d y ,d z }). 




Figure 2: Plot of the pseudo-parametrization used, with its 
distortion through space [Eq. ( |10[ ) with a = 1, k = 2 and 
6 = 5] 



Using these to transform the Gross-Pitaevskii equation 
[Eq. 0] into the local coordinates around a point on the 
vortex line r , assuming that ft is in z, and assuming 
that the vortex line has no velocity in z, the equation for 
the behaviour near the vortex core becomes 
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dx 
V* = ih 



dy 
d^ 

~dt 



V tr (r ) 



(17) 



where the trap terms have been expanded to 1st order 
(order r c ), r is the position vector in the local coordinate 
system, fl is the rotation vector and V is the velocity of 
the vortex line at Tq. 
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This equation describes the behaviour of the conden- 
sate near the vortex core, in terms of the properties of 
the vortex core itself. Hence, solving this equation gives 
the desired inner region behaviour of the BEC. 



IV. THE BEHAVIOUR NEAR THE VORTEX 
CORE 



The steady state solutions of Eq. ( 17) describes a con- 



densate containing a vortex with standing waves along 
the line of the vortex. Hence to determine a relation for 
these standing waves, the inner equation to solve is 
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2M \ dx dx dz dy dz 

+V tr (r ) * + g|*| 2 * 
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a 9_ +/3 _d_ 
dx dy 
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V tr (r ) 



(18) 



where V = V + ($7 x r ). 

The small nature of the wave perturbations give four 
small parameters (assumed to be of similar order) to use 
for a perturbation expansion: k, a, /3 and r c . The Oth 
order terms, of this expansion describe a condensate with 
an unperturbed vortex: 
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The 1st order terms account for the perturbations on the 
vortex line: 



y/\JI(ffj — Vtr(v )]/g. This simplifies Eq. ( 20 ) to 

h 2 
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where d z i = d z + ad x + (3d y , and |VVt r (ro)|i and V- 
are the trap gradient and vortex line velocity in the i 
direction. 

Assuming that the wave function has the form 



* = [|* | + D(p, z) + X (p, z) cos{4>) 

+ C(p, z) sin(0)] e i i't>+^{p,z) cos(</>)+zA(p, 2 ) S in(<A) ^ 



(22) 



makes * = 1*0 1 ^ and 



*i = [D(p, z) + x(p, z) co$((f)) 

+ C{p, z) sin(0) + i \9 \v{p> z ) cos(0) 
+ t|*o|A(p,z)sin(^)]e < «* l (23) 

where q is the winding number of the vortex. 



This reduces Eq. (21) into five coupled differential 



equations describing the behaviour of the perturbation 
functions D, x, 77, and A: 
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V tr (r ). 
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Only the far from core behaviour of the inner so- 
lution is used for matching. Hence, the large p = 



\/ x 2 + y 2 behaviour of Eq. (20) is required. In this 



regime |^ol should be cylindrically symmetric, and can 
be approximated by the Thomas-Fermi solution I^tfI ~ 



Id, <9L> d 2 D q 2 D 4Mg\^ \ 2 D 
(p^-)+ ' 



p dp dp dz 2 p 2 



h 2 

2Af|*o| z 

h 2 



d z ,V tr (r ), (24) 



1 _L( + _ (g 2 + l )x _ 2g|*o|A 

p dp dp dz 2 p 2 p 2 

4Mg\^ Q \ 2 x _ 2M|tf ||W tP | a p 



h 2 h 2 

2Mq\^ Q \V; d|* | 



hp 



dp 



2a 



d 2 |*o| 
dp dz 



(25) 



5 



I A( dc ) i 92 c (g 2 + 1 )c 
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h dp p dz 



(26) 



(27) 



order trap radius R±; therefore the change in \^>o\ 
with respect to p (varies on r c ) is small compared 
with IV'ol/Pi which is approximately constant, and 
so has little effect. 

• Kinetic energy terms, from the waves density, can 
be ignored. In the Thomas-Fermi regime the ki- 
netic energy from the waves density is assumed 
negligible, giving the well known inverted parabola 
structure (c^tf! 2 = PTF — Vt r )- As this behaviour 
is expected in the matching region, the Laplacians 
(V 2 ) of Eqs. ( |24| ), p5| and ([26]) are approximately 
0, removing the differential behaviour from these 
equations. 



Applying these assumptions and rearranging Eqs. (24 1, 
p5l), (p6l), d27|, and p8l), we find 
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dz dz p 2 
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dz 
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These equations separate into two connected groups: 
equations describing perturbations on condensate den- 
sity [Eqs. (24), p5| ), and ([26])] , and equations describing 
perturbations on the condensate's phase [Eqs. (271 and 
(28)]. Different groups having different assumptions and 



transformations that simplify the set of equations. 

• Substituting rj and A with fj + pV^M/H and A + 
pV^M/h, removes V dependence from the describ- 
ing equations. As p is a solution to d 2 f(p) + 

dpf(p)/p~f(p)/p 2 = 0: an arbitrary constant mul- 
tiplied by p can be added to f(p) without changing 
its final solution; consequently, this constant should 
be chosen to best simplify the equations, removing 
the V dependence. 

• Derivatives of the phase terms with respect to z 
can be omitted (and by association fj and A). The 
gradient of the phase defines the condensates ve- 
locity, variation in z therefore inducing flow in z. 
This kind of flow is caused by either the addi- 
tion of sources and sinks or vortex rings to the 
BEC. This derivation only considers the behaviour 
of a 'straight' vortex in a condensate, and therefore 
should have no velocity in z. 

• The phase terms d p \^fo\d p r} and d p \^>o\d p X can be 
removed. In Eqs. (27) and (28) there are two terms 



involving the 1st derivatives of the phase with re- 
spect to p, having coefficients \tpo\/p and 9 p |^o|- As 
mentioned previously, at large p, \&o\ 1S approxi- 
mately I^tfI, and varies on length scales of the 
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2qC 

r.2 



dp 2 
2qad\^ 
p dz 



A 
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(29) 



(30) 



(31) 



(32) 



(33) 



Equation (29) solves for the perturbation D(p,z), be- 
ing the first order correction to the local Thomas-Fermi 
profile, i.e. correcting for the changing condensate den- 
sity in the z direction. 

The remaining density perturbations, Eqs. (30) and 
(31), are linear. These equations contain appropri- 
ate Thomas-Fermi, curvature and i deviation correction 
terms and a term that depends on a phase perturbation. 
As a consequence, these will be solved once the phase 
perturbations are known. Using these to remove x and 



C dependence from Eqs. (32) and (33), two independent 
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differential equations for T) and A are found: 
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Comparing 

fi 2 g« gtjgoj 



the relative 

?i 2 ga 3 2 |W | 



terms, 
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Mg|* |p 2 ' 
K 2 g/ 3 d 2 [gg I 



2M s |* "Fp 2 dp ' Mffhf^fp 2 <9p dz ' A/g|* |V 9p dz 

are much smaller than their corresponding counterparts 
by at least an order of 1/R± (especially as p gets large). 
Consequently, these terms do not significantly affect the 
far from core behaviour and can be omitted, turning 



Eqs. (34) and (35) into 
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q\Wtr\ x 

g |*TF f P 



(37) 



where |\&o| has been replaced with |*tf|- 

Given that |\I/tf|, |VVt r | and are constant in the 
inner solution expansion, these equations are Euler ho- 
mogeneous equations, and have solutions of the form p n . 
Solving these equations, and making them non-divergent 
at p — the solutions become 



<z|w tI 



qP dV tr 



2g\V 



TF 



TF I 



dz 



pln(p) 



M TTI 

+ -j-pV* + pA 




where A and B are unknown constants that can depend 
on \V TF \, |W tr |, k and ^f. 

Therefore the pha°" * lw> ; • '•><!<••• >• 
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(z x 



' dz 
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~ z x VVt r (r ) z x t dV tr 
g\*TF\ 2 ' g\^T F \ 2 dz 



,r 



r=r 



(41) 



as t = {a, ft, 1}, and E = {A, 5,0}. The dependant 
terms in E are written down here to emphasise that the 
constant can depend on |*tf|, |VV tr |, k and with 
the structure of the dependants chosen because of the 
structure of the other terms in Eq. (41 ). 



Using Eqs. ( 38 ) and ( 39 1 , the density perturbations are 
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2g\q T F 
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(44) 



These density perturbations are independent of the vor- 
tex core's velocity (V), being, to lowest order, corrections 
to the local Thomas-Fermi profile. Hence these solutions 



7 



already match the Thomas-Fermi density profile of the 
outer region. This indicates that the perturbations on 
a vortex line predominantly affects the inner phase of 
the condensate [Eq. (41 (], which through matching to 
the outer regions phase, should give the desired vortex 
relation. 



V. THE PHASE OF THE OUTER REGION 

Far from the vortex core, the BEC's wave function 
density has the Thomas-Fermi profile, and its phase (S) 
solves 

V • (\V TF \ 2 VS) - • (r x V) \^>tf\ 2 = 0, (45) 
v ' h 



z-Vx (VS)=2irqSW(p-p ), 



(46) 



where S^(p — po) is the delta function in p and <f>. 

These equations are linear, and therefore allows S to 
be solved in two separate components: one dealing with 
rotation [the particular solution of Eq. (45), Sq] and one 
with everything else (S c ): 

V ■ (\V TF | 2 V5o) - -z-fl ■ (r x V) |* T f| 2 - 0, (47) 



V x (VS n ) = 0, 



V- (|*tf| 2 V5 c ) =0, 
i ■ V x (V5 C ) = 2nqS<- 2) (p - p Q ). 



(48) 



(49) 



(50) 



Assuming fl = flz Eqs. (471 and (48) have the solution 

,2\ 



5* 



M{<4 



h (w2+ w 2) 



flxy, 



(51) 



where oj x and u> y are the trapping frequencies in x and 
y, respectively. 

Transforming into the local coordinate system of the 
line (r — > r + r ), Eq. (51 1 becomes 



Mf , , o VVtr(r ) x fl \ 

— (0 x r ) + 2 ; ' I • r. 



(52) 



S c is not as simple to solve, with no exact solution ex- 
isting for arbitrary p a . The gradient of S c , however, can 
be determined and consequently compared with its in- 
ner solution equivalent. To do this VS C is written as 
-g(V x (|'I' tf | 2 <I>z))/|<I't F | 2 , so that it automatically 



satisfies Eq. ( |49| ) , and turns Eq. ( 50 ) into 
i ■ V x (V£ c ) 



V • (VS C x 2) 
V- (V_L$ + $V_Lln(|* TF | 2 )) 

2^ 2 \p-p ), (53) 



which through rearrangement becomes 



27r<5( 2 )(p- p ). 



(54) 



Transforming Eq. |54j) into the local coordinates of the 
vortex line, and taking V^e - ln to be approximately 
e -l n |*TiHy2 V tr /2g\^ TF \ 2 , it becomes 



TF | 



2^(p)e ln l* r - 



'I 



(55) 



In this outer region there should be no flow in the p 
or z, and therefore $ is purely a function of p. Further- 
more, J. R. Anglin [2g a showed, in a complete hydrody- 
namic calculation, flow only arises from the vortex's vor- 
ticity, thereby indicating that the homogeneous solution 
to Eq. (55) is irrelevant. Hence the solution to Eq. (55) 




(56) 



where K is a modified Bessel function of the second kind 
and Vtr and gl^TiH 2 are implicitly functions of ro- 
For small p, 



(57) 



where c 



0.577... is the Euler constant, and 



Vj_Vt r /2g\^TF\ 2 has been approximated by 1/R± (as 



V 2 ± V tr « Mw} and g\^ T F\ 2 



Using this 



I* 



TF I 



V x |* TF | 2 ln 



e p 
2~R~\ 



p 



<vV x ln ( ) z 
gin V x 1h|*tfI~z 



2 s 



<?|<J< TF | 



z x VV tr (r ), 



(58) 



having used the identity V x A(r)z = — z x Vj4(r). 
Therefore, the outer phase is fully described by the 



Eqs. (58) and (52). Through matching these equations 



with the inner phase [Eq. (41 )] a relationship for the ve- 



locity of the vortex line is determined. 



VI. MATCHING THE INNER AND OUTER 
SOLUTIONS 

To match the outer and inner solutions, the inner phase 
must be divided into the same Sn and S c structure as the 



outer solution. Assuming that the core velocity V can 
be written as a rotational and other V c component, 



Eq. (41) separates into 
M 



Sa 



(V n + nx T ) ■ r, 



and 

S r , = 



h 2 



2g\* 



TF 



z x VF tr (r ) + (z x t) 



r E 



Kb 



z x VV tr (r ) i x t 9V ir 



TF| 



.9 1* 



TF I 



92 



• r 



(60) 



which can be compared to the outer Sq [Eq. ( 52 )] and S c 
[Eq. (58)] to find the matching condition for the velocity 



of the line. 



Comparing Eq. (|59| with Eq. (|52|), shows that Vq is 

(61) 



W tr (r ) x SI 



To determine V c , however, V5 C of the two regions 
must be matched. Using Eq. (60), V5 Cj for large p is 



-0 + T V c r - ^— • In p + p ® p] 

p h 2 



q zx W tr (r ) 



2g\V TF \ 
q (z x t) 



[ln(p) + p ® p] 



t\ dVtr\ 



2g\q TF \ 



E, 



(62) 



as for large p, - p and Vr = I (where I is the identity 
matrix) . 



Matching Eqs. (62) and (58), at the vortex core radius 



r c , the equation for V c becomes 
hq nb 



V c (r ) 



[ln(r c ) + p ® p] 

pn(r c ) + p ® p] 



2M 

S?zx VV tr (r ) 



+ 



2A/ 5 |* TF r 

^(^xt)^| r=ro 



2M 5 |<J< TF | 



[ln(r c ) + p <g> p] 



Combining V c [Eq. d63b] and V n [Eq. (fSlb] the rela- 



tion for velocity of the vortex line is 



(59) v(r ) = 



[ln(r c ) + p(g) p] 

[ln(r c ) + p <g) p] 



2M 

frg z x VV tr .(r ) 
2A/g|* TF | 2 



2M 5 |* rF | 



[ln(r c ) + p <g> p] 



E' 



Kb 



z x VV tr (r ) z x t 9Vt, 



9 Y&TF 



5 I* 



dz 



,r 



TF 



[z x W ir (r )] +2 



W fr (r ) x n 
(64) 



where p is effectively the p vector in the condensates co- 
ordinates (the difference between p of the vortex line and 
p of the condensate is negligible at large distances). 

This equation indicates how the local velocity of the 
vortex line behaves in response to the lines position and 
curvature, and therefore reveals the local behaviour of a 
vortex line. 



Equation (64) also transforms to give a system of dif- 



ferential equations that describe the complete vortex 
behaviour. As the vortex line is a one to one func- 
tion in z, the vortex line's shape and motion are fully 
parametrized through two functions, p(t,z) and cf>(t,z). 
This allows any point on the line to be represented by 
r = p(t,z)p + <p(t,z)4> + zz, the velocity and tangent 
vectors to be 



Velocity = V 
tangent vector = t 



dp , d(j) 



dt 



p + p 



dt 



dp , d<f) 



dz 



P + p—4> + 

dz 



(65) 



(66) 



and Kb (« i x k, where k = d 2 Vo) to equal 



Kb 



d 2 <l> +2 ^P^_ 
dz 2 dz dz ' 



E' 



- z X VVt r (r ) i x t dV tr 



where E' 



,g|* TF | 2 '<?|*tf| 2 9 
hq\n(^A 

[zxVy ir (r )], 

Mg\W T F\ 2 
-KE/M. 



,r 



cPp 
dz 2 



dip 



(67) 



(63) 



where p and <fr are implicitly functions of t and z, and z 
has been approximated to be the arclcngth of the curve. 

As a result, the differential equations describing the 
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full vortex's shape and motion are 
hq [ln(r c ) + 1 



dp 

at 



2M 
Hq [ln(r c ) 
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2Mg\V TF \ 



dz 2 

ldV tr 



dz dz 
d4> 8V tr 



2tt dV tr 

P V\Vtr ~W 

p-FJ, 



p dcj) dz dz 
pMg\^TF 



dVtr 

2 dcj) 



(68) 



d(f> ftqln(r c ) 
'7ft = 



2M 



(Pp 
dz 2 



P 



dz 



+ 



ftgln(r c ) ( dV tr dpdV t7 



2Mg\^ TF \ 2 V 9 P 

2Vi dv tr h i ln 



v 2 ± v tr 



dz dz 

dV tr 



c r 
2R_ 



dp Mg\^TF\ 2 dp 



(69) 



The solutions of Eqs. ( 68 I and ( 69 ) describes a vortex's 



shape and motion, and shows the steady state structures 
and waves the vortex supports. However, before these 
structures are determined the unknown constant E' needs 
to be found. 



VII. DETERMINING E' 



Before investigating Eqs. (68) and (69) E' needs to be 



determined. This can be done through comparing specific 
results from these equations to a previously determined 
physical scenarios. 

Assuming that the trapping and curving dependence 
in E' can be treated separately, two different scenarios 
need to be considered: one indicting a straight vortex's 
reaction to the trap confinement, and one describing be- 
haviour of wave perturbations on the line in the absence 
of confinement. 

In a pancake shaped condensate it has been shown 
that the vortex lines are straight, while for more cigar 
shaped traps the line tends to bend at the edges [15] , 
Such straight lines are easy to work with mathematically 
and so the energetics and behaviour of these lines have 
been investigated [12 E] . These investigations revealed 
that a off centred straight line vortex precesses around 
a cylindrical condensate in the direction of the vortex's 
rotation with a precession frequency 



3ft In 



dt 2MR\ (1-4 



n. 



(70) 



Applying similar conditions to Eqs. ( 68 1 and ( 69 ) (dis- 
carding z derivatives and setting Vt r = M(Jjjr/2 + 



Muj 2 z 2 /2) and setting g\^ TF \ 2 « \Mu)\R\ (l - jjr) 
the describing equations become 



dp 

dt 



dt 



H In 



0. 



2 p( p 



_ #x(i-£) 
(«f)> 



,0,r 



(71) 



MRj_ 1 



Qp 



E' 



0. 



2p<f> 



,0,r 



(72) 



When matched to Eq. (|70|) these equations become 
p-E' 



0. 



2p4> 



,0,r 



0, 



E' 



0. 



2p<t> 



41-1 



,0,r 



ft In 



mrUi-£ 



(73) 



(74) 



This suggests that the trap dependence in the constant 
term has the form 

, / c Z X WV tr Z X t dV tr \ 

E nb, s-, ^— — ,r 



g\y TF \ 2 ' g\^TF\ 2 dz 

^ ln (4^fe)( 2xV ^) 



2Mg\y TF \ 



+ F Kb, r , 



(75) 



where F is an unknown constant that describes the cur- 
vature behaviour. 

Note that this condition, though intending to contain 
all trap dependence, does not contain information on 
the lines reaction to the z trapping of the condensate. 
This indicates that another condition relating specifically 
to this term is necessary (visible through the In term). 
Unfortunately a suitable comparison was not found and 
therefore this term could not be calibrated. 

To determine the curvature related behaviour of the 
constant, the behaviour of a perturbed vortex (in a un- 
trapped condensate) needs to be known. It has been 
shown, theoretically, that vortices in BECs support he- 
lical wave perturbations [TU] and therefore are an ideal 
example. These waves rotate in the opposite direction to 
the rotation of the vortex, and in a uniform condensate 
has a dispersion relation (for large wavelengths) of 



ftfc 2 
2M 



hi 



1 



fc|< 



(76) 



where w is the wave frequency, k is the wave number 
and r c is the vortex core radius. Hence by comparing 
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this dispersion relation to the equivalent derived from 



Eqs. (68) and (69) the form of F can be determined. 



A helical wave generally has the form x — a sin(/cz — 
cut + (f>o) and y = acos(kz — cut + </>o), which when 
transformed into cylindrical coordinates is p = a and 
tf> = kz — cut + 4> (a being the wave amplitude and <f>o 
being an arbitrary phase constant). Substituting this 
parametrization into Eqs. (68) and (169]) , they become 



= _ fi[ln(r e ) + l] / fc2 ^ 

2Mg\y TF \ V - 



(77) 



and 



hak 2 ln(r c 
2M 



-ak .(f), r 



(78) 



Hence, the conditions F satisfies are 



p-F(-ak 2 (f),r) = 0, 
F(-ak 2 (f),r)=a—H\k\r 2 c ) 



(79) 
(80) 



Since all the wavenumber dependence in F must come 
from the Kb term a suitable form becomes 



tion VI I , this equation gives 

hqp [ln(r c ) + 1] 90 dV tr 



dp 

at 



2Mg\f TF \ 2 dz dz 



hq 



In 



2M 



d 2 (f> dpd(f> 



dz 2 



3hq 



In 



(*) 



dV t , 



and 



d(f) 

'at 



ApMg \^tf\ 

n av tr 

J pV 2 ± V tr d<$> ' 



hq\n(r c ) dpdV tr 
2Mg\^ TF \ 2 dz dz 

hq\n ' 1 



dz dz 



(83) 



3%lnjj£) dVtr 

4Mg\^ TF \ 2 dp 



pid^Y-dlp 



O 



2M 

dVtr 



(pp 

dz 2 



V\V tr dp 



(84) 



From these non-linear equations the complete be- 
haviour of the vortex line can be determined, however, 
this paper will only look at some general and simple ex- 
amples of these equations. 



hq In 



2M 



Kb. 



(81) 



Substituting the full form of E' into Eq. ( 64 1 the ve- 
locity condition becomes 



V(r ) = 



ln(r c ) + p <g) p] 




Hq(z x VV fr (ro)) 
2Mg\V TF \ 2 

2 vVtr(r ) x n 

V 2 ± V tr (r ) 



P 

(82) 



This is an almost complete relation for the velocity of 
the vortex line, missing the calibration of the d z Vt r term 
[as evidenced by the units in the multiplying logarithm 
ln(r c )]. As a suitable analytic form for this comparison 
could not be found, and it does not affect any of the latter 
calculations performed, this issue will be left for future 
work. 

When transformed into differential form (as per Sec- 



VIII. WAVE INVESTIGATIONS 

A. Qualitative comments 

These equations give deterministic properties of a 
slightly perturbed vortex in a BEC. However before look- 
ing at a special case of this, it is worth looking at quali- 
tative behaviour from Eq. ( 82 ) . 



The vortex contains angular momentum in z. There- 
fore, as suggested by A. Fetter [7j, the vortex can be con- 
sidered as a gyroscope spinning in z, responding to any 
force by moving perpendicular to it. In Eq. (82) there 
are 3 such forces: a trap gradient force (— VV^TTro)), a 
curvature force (m, b ps i x n), and a z-trap deviation 
force (d z Vt r (ro)t). Each force contributing to the rela- 
tions final form. 

Similarly, the rotation of the condensate f2, also creates 
gyroscopic motion. This rotation causes any point in the 
trap to move perpendicular to the trap gradient force, 
and therefore also pushes the vortex line accordingly. 

Equation ( 82 ) also contains a collection of terms that 
'amplifies' the velocity of the line in the p direction. 
These terms are additional restoring forces, which at- 
tempt to return the system to the the centred straight 
line vortex state. 



Hence the qualitative behaviour depicted by Eq.(82) 



gives understandable insight into how the vortex line 
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moves; therefore, indicating that it is practical to use 
it to quantitatively determine the structure and motion 
of the vortex line in a condensate. 



B. More realistic helical wave behaviour 

As previously mentioned, helical waves can exist on 
vortex lines in a BEC. Hence it is relevant to consider 
the behaviour of such waves in the presence of a trapping 
potential. Assuming a cylindrical trap (Vt r = Muj 2 _p 2 /2, 
where lo z — 0) and that the helix wave takes the same 
formed used in Section VII (p = a and <j) = kz — u}t + 4>o), 
Eq. ( 83 ) is solved and Eq. ( 84 ) becomes 



3fi?ln(^) 



hqk 2 In 



l 

r e |fc| 



2MR\ 



1 - 



2M 



n. (85) 



This equation shows that the precession helical wave in 
a cylindrical trap is a combination of the precession of a 
straight line vortex and the motion of a helical wave in 
an uniform condensate. 

In a numerical simulation by T. Simula et al. \29\ 
[50] . it was shown that the dispersion relation of large 
wavelength helical waves in a cigar shaped trap could be 
represented as 



uj 



^ fc2ln (^) 



2M 



SI, 



(86) 



where uj is a constant precession frequency. The similar- 
ity of these two equations is striking, and upon plotting 
(Fig. |3| , the similarities are still visible (in the change to 
to to the changing particle number), however the specific 
values of the constant vary by about ±0.1(u;/wj_) to the 
numerical value. 

This deviation is not surprising. Equation ( |85[ )'s con- 
stant precession is that of a purely straight line vortex. 
This suggests that Eq. ( 85 ) is representative of pancake 



shaped traps where such vortices occur. T. Simula's work 
however is in a cigar shaped trap and therefore the vortex 
line, without wave perturbations, tends to bend, poten- 
tially changing its precession. Therefore to accurately 
depict this behaviour the trapping of the condensate in 
z, and how it distorts the vortex, needs to be considered. 



IX. PREVIOUS SIMILAR WORK 

A similar derivation has been carried out by A. 
Svidzinsky and A. Fetter [25 a with the aim of determining 
the normal waves that a vortex line can support. This 
work performed a similar matched asymptotic expansion 
to determine a relation for the vortex lines velocity. 

During this derivation, an extra term appeared in the 
outer equations that was related to the transformation 
into the local coordinates of the vortex line, the specifics 




Figure 3: Plots of the dispersion relation, Eq. ( |85| , for dif- 
ferent particle number (N) for a condensate of Rb atoms, 
with trapping frequencies uj± — 2n x 11.8 Hz, and uj z = 
2tt x 98.5 Hz, and vortex core radius r c — 0.13pm. 



of which were not explained in the paper. This gave an 
outer solution of 



-e 2 K [ p< 



V 2 ± V tr 




(87) 




2g\V TF \ 
which was approximated to be 



for small p. This form is almost identical to that derived 
in this paper, having effectively ignored the additional 
exponential in their outer solution. 

Then to compare the outer and inner solutions a 
pseudo-$ was created for the inner phase solution. This 
pseudo- ( I > was linked with the inner solution phase S 
through the same definition of <!> used in Section [Vj but 
ignored some terms, of comparable order to those kept, 



1 



|VxVtr| 

2 \g\*TF\ 2 



MVn 

MEp) + — —x. (89) 
nq 



It is worth noting that the largest term of this inner $ 
is xln(p), while the largest term of the outer $ is ln(p). 
This indicates an inconsistency in the expansion of the 
outer solution (the e^ in Eq. (87) allowing the outer to 
also expand to 2;ln(/?)). 

Through matching these <E>s the velocity condition 
found was 



V(r ) 



qhi x VV tr 
2M 5 |* TF | 2 

qhuh 



In 



In 



2M 

VV tr X ft 



Wc 







K 2 










■t) 











(90) 



This condition gives similar wave results within certain 
limits. Specifically both equations replicate the long 
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wavelength helical wave dispersion relation but generally 
give different results and dispersion relations. 

Furthermore, in their derivation the effect of t devi- 
ating from z was not considered, but later implemented 
through the t x VVt r term. This procedure gives the 
same d z Vt r structure seen in Eq. (82) however also adds 



an additional artefact in the i direction of the velocity 
equation. 



X. CONCLUSIONS AND FUTURE 
DIRECTIONS 

Through this derivation, a robust relation for the mo- 
tion and shape of a vortex line in a BEC was determined 
[Eq. (82)]. This relation was found by determining the 
behaviour of the condensate near and far from the vor- 
tex core separately, and ensuring that these two solutions 
matched in the overlapping region. 

This relation was shown to describe well the local vor- 



tex's behaviour, and could be transformed into a set of 
differential equation that describes the vortex's overall 



structure and motion [Eqs. (83) and (84)]. Using a sim- 
ple calculation, these equations reflected wave behaviour 
seen numerically. 

This paper only briefly explored the behaviour of 
Eqs. (83) and (|84j) , and so further investigation needs 



to be performed, looking at different vortex line struc- 
tures and trapping geometries. However before this can 
be done an analytic description of how the trapping in 
z affects the line's behaviour is needed to calibrate the 



final term in Eq. ( 82 ) 
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